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Abstract 

In this work the implicit function theorem is used for searching local 
symbolic resolution of differential equations. General results of existence 
for first order equations are proven and some examples, one relative to 
cavitation in a fluid, are developed. These examples seem to show that 
local approximation of non linear differential equations can give useful 
informations about symbolic form of possible solutions, and in the case 
a global solution is known, locally the accuracy of approximation can be 
good. 

1 Introduction 

The use of Implicit Function or Dint's Theorem for local resolution of differen- 
tial equations is not new (see [7] for an historical example of application due 
to Bernoulli, and [9] for an application of theoretical existence for a non lin- 
ear system of PDEs). But in the technical literature this method is very little 
widespread and limited to some specific situations. Moreover, for some physical 
and engineering aims, often usual numerical resolution of differential equations 
is not very useful, because the knowledge of an analytical dependence from 
physical or geometrical parameters is required (see e.g 3 ) for the problem of 
searching analytical dependence for resistence coefficients in flows). 
In this work some statements about local symbolic resolution of ordinary dif- 
ferential equations are established, and some examples, two very simple but 
illustrative the technique, and other two more technical and useful, are devel- 
oped. 
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2 ODEs and Dini's Theorem 



Let [a, b], a < b, an interval of R. If we denote with y = y(x) a C n ([a, &]) 
function, for some n integer, n > 0, defined on [a, b] and with j/' -' = y^ k \x) 
its A;-th derivative, an ordinary differential equation (ODE) is an expression like 
this: 

F(x,y(x),y^(x),...,y^) = Vxe[a,6] (1) 

where F is a function F : [a, b] x R™ +1 — > R. Let d x F the partial derivative 

of F respect the x variable in [a,b] and duF, < k < n, the partial derivative 

respect a suitable variable of R n+1 domain. 

Suppose now F E C m ([a, b] x R n+1 ) for some integer m > 0. 

If T = (xg,pg,...,p°) € [a, 6] x R™ +1 is a point such that F(T ) = 0, and 

dkF(To) for a particular value of k, then from Dim's or Implicit Function 

Theorem there are a constant 5 > 0, a costant M > 0, a neighbour Uq of 

(pp, ■■■,Pn) and a unique function : [a, 6] x 1" — > R, S C m ([a, 6]), such that 

F(x,p , -,Pk-i,<t>(x,p, ...,p k -i,Pk+i, ...,p n ),Pk+i, -,p n ) = (2) 
V(x,p , ...,Pfc-i,Pk+i, —,Pn) € (x - <5, x + <5) x J7 

0(a;8,p8,...,pg_ 1 ,pg +1 ,... J p°)=pg (3) 

|0(^O,PO,-,Pfc-l,Pfc+l,-,Pn) <M (4) 

and 

90 = _^T 

dp, 9 fc T 1 ' 

For a complete proof, see e.g. 0|. 

Let we suppose n = 1, so the differential equation is F(x, y(x),y'(x)) = 0. If the 
equation has order greater then 1, it is possible to write it as a system of order 
1 differential equations (see e.g. [Sj), and Dini Theorem has a natural extension 
to systems of equations. 

Theorem 1 Let F(x, y(x), y'(x)) = an ODE 7 with F e C 1 ([a, 6] x R 2 ) a real 
function F(x,p, q), and let (xo,Po, qo) such that xq € (a, b), F(xo,po, qo) = and 
d q F(xo,po, qo) ^ 0. Then exist S > 0, p > 0, a function u £ C ([xo — 5, xo + 5]) 
and a function <fi G C°([xo — 5, xo + 5] x [po — /O,po + p]) swcft i/iai 0(xo,po) = <Zo 
and 

n'(x) = c/f>(x, w(x)), F(x, ti(x), 0(x, u(x))) = Vx € [xo — S, xq + S] (6) 
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Proof. From Dini Theorem exist Si > 0, p\ > and a function G C°([xo — 
Si,xo + Si] x [p — Pi,Po + Pi}) such that 

F(x,p,<f>(x,p)) = V(x,p) € [x - Si,x + Si] x [p - pi,p + p x ] (7) 

We can consider the differential problem 



Being <f> a continuous function, from the theorem of Peano about existence of 
solutions for first order ordinary differential problems (or by local existence the- 
orems) (see 0,111), exist 82 with < 82 < 8\ and a function it G C 1 ([xo — 82, xo + 
62]) (local) solution of (JSJ. Being u a continuous function and u(xq) = Po, exists 
a < S3 < §2 such that u{x) G [po — PiiPa + Pi] Vx G [x — <5 3 , x + <5s]). 
Then, from Q, if x G [xq — 83, xq + S3}) C [xq — ^i,xo + <5i]), we have 
F(x, u(x), <f>(x, u(x))) — 0, so with the choice 8 — 83 and p = p\ the proof 
is completed. □ 

Therefore, under suitable hypothesis, ODE F(x, y(x), y'(x)) = has at least 
local solutions. But the theorem says too that the original ODE can be locally 
written in a new form, i.e. in the expression y'(x) = <f>(x,y{x)), that is the 
normal form. In this sense, the term locally means that there is a neighbour of 
(xo,PO) 9o) where the original ODE can be written in the new form, and a local 
solution u exists such that u(xq) = po and u'(xq) = qa- 

Example 1. Consider the ODE 



in the interval [a, b] where a < and b > 0. The function F is F(x,p, q) = 2p — q, 
with F G C°° ([a, b] x M 2 ) and d q F = -1. Let x = and p = 1; then 
qo = 2. From theorem 1 exist two functions u and <j> such that </>(0, 1) = 2 
and 2p — (f>(x,p) = in a neighbour of (0, 1). Then (j)(x,p) = 2p, so from 
u'(x) = 2u(x), that is it(x) = ce 2x . The constant c is determined from the 
condition u(0) = p a = 1, hence c = 1. 

In the case d q F(xo,po,qo) = 0, we can consider the possibility of a local 
expression like y(x) = ip(x,y'(x)), as stated in the following 

Theorem 2 Let F(x, y(x), y'(x)) =0an ODE, with F e C 1 ([a, b] x R 2 ) a reaZ 
function F(x,p,q), and let (xo,po,9o) such that xq g (a,b), F(xo,po, <7o) = 
aKc? d p F(xo,po, qo) 7^ 0. TTien exist (5 > 0, p > and a function ip G C°([xo — 
J, xo + S] x [go — p, qo + p]) such that ip{ x o, Qo) — Po an d 

F(x, ip(x,q),q) = V(x, g) G [x - 8, x Q + 5} x [q - p, q + p] (10) 

If exists a solution u G C ([xo — a, xq + a], with < a < 8, for the problem 



y'(x) = <f>(x,y(x)), y(x )=Po 



(8) 



2y(x) - y'(x) = 



(9) 



y(x) = ip(x,y'(x)) 
y'(x Q ) = g 




(11) 
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then 



F(x, u(x), u (x)) — Vx G [xq — a, xq + cr] 



(12) 



Proof. The first statement is a direct consequence of Dini's Theorem. Let then 
u a solution of the problem (|llfl . with o\ < S. From the fact that u' is con- 
tinuous in the interval [xo — a±, xq + ci], we can choice the constant a so that 
u'{x) £ [go — P, 1o + j°] for every x £ [xq — er, + cr]. Then, from (|10fl follows 
that F(a;, u(x), u'[x)) = for every x £ [xq — a, xq + a]. □ 

Therefore, a local solution for the new differential equation y(x) = ij){x,y' (x)) 
is a local solution for the original ODE F(x, y(x), y'(x)) = 0. 



Example 2. Consider the non linear ODE 

y'(x) 2 +y(x) 2 



1 = 



(13) 



The function F is F(x,p, q) = q 2 + p 2 - 1, with d p F = 2p and d q F = 2q. If 
we consider the point T = (|, 1,0), we have d q F(T ) = and d p F(T ) = 2, 
so from the first statement of theorem 2, if \q\ < 1 exists the function ip — 
ip(x, q) — \/\ — q 2 such that (locally) p = \Jl — q 2 . Then we can consider the 
problem 

y/l - y'{x) 2 Vx £ [0, tt] 



y{x) 
y'(f) = o 



(14) 



The function u(x) — sin(x) is a solution of this problem, with |w'(x)| < 1 as 
required, and is a local (but in this case global too) solution of the original ODE. 
Note that finding the solution of the original equation implies the use of complex 
variables: the internal command DSolve [y' [x] ~2+y [x] ~2-l==0,y [x] ,x] of 
the software Mathematica gives the integrals 



u(x) = 1 — 2sin 2 



~(x-ic) 



v(x) = 1 — 2sin 2 



\{~x-ic) 



(15) 



where i is the imaginary unit and c a complex constant. We obtain our previuos 

j|7r and c = 



local solution respectively for, e.g., c — i\~K and c = — i§7r. 



3 Local symbolic approximation for an ODE 

From implicit function theorem, if F has continuous partial derivatives, using 
Taylor expansion near the point {xq,pq) (or (xq, qo)) and the chain rule for the 
derivatives of composite function it is possibile to write the following expression 
for the function (or ip): 



(x,p) = (j)(x ,p Q ) + [d x (j)(x ,Po)] (x-x ) + [d p cj)(xo,Po)] (p~Po) 



(16) 



+ 



-dl x (j>{xo,po) 



(x - x ) 2 + 



1 . 

2°pp 



dlMx ,p ) 



(P - Paf + 



with, e.g. for ^-derivation, 



Q x p 

d x (j>{x ,p Q ) = --^-=(x ,Pa,qo) (17) 

d xx <f>{x ,Po) = (18) 

where the derivatives of <j> are evaluated at (xq,pq) an d those of F at (xq,pq, qo)- 
The expansion gives a local analytical approximation of the original ODE. If for 
this equation a global solution doesn't exist or it is not given by elementary 
functions, the expansion could be a more simple equation to solve. From pre- 
vious theorems, a solution of this equation is an approximate local solution of 
the original ODE. 

If F S CP-fab] xl 2 ), F(x ,p ,q ) = and d q F(x ,p , q ) ^ 0, from 
previous identities the first order differential equation F(x,y(x),y'(x)) has the 
following local approximation of first order. 



y'(x) =q + 



d x F 
d q F 



(x ,p ,q ) 



•< - {) 



(x ,Po,qo) 



d x F 



L d q F 



{x ,p ,qo} 



d p F 
d q F 
d P F 

g-p{xo,Po,qo) 



Po - 

y[x) 



(19) 



If d p F(xo,Po, Qo) = 0, the solution of this first order ODE is 



Qo 



d x F 
d q F 



xo 



ld x F 



2 d q F 



otherwise its general solution is (see e.g [[]) 



y(x) = Pa 



d x Fd q F d x Fxo + d q Fgo 8 X F 
[d p F] 2 d p F 8 P F 



x + c e 



(20) 



(21) 



where c is a constant and the partial derivatives are all calculated at point 
(xo,po,qo). Imposing the condition y(xo) = po we obtain 



[d x F + d p Fq ] 3 q F e 



0„F X0 



[d P FY 

while, with this value of the constant c, we have y'(xo) = qo as required. 
So we have proven the following 



(22) 



Proposition 1 In the hypothesis of Theorem 1, the ODE F(x,y(x),y'(x)) = 
has a local solution, approximated to first order, of the form 



y{x) 



a + b x + c x 



if d p F(x ,p ,q ) = 



a + b x + c e d ° x ifd p F(x ,p o ,q a )^0 



(23) 
(24) 
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where ao, bo, cq and do are constants depending from (xo,po,qo)- 



Example 1 bis. For the ODE 2y(x) — y'[x) — 0, being F(x,p,q) — 2p — q, at 
point (0, 1, 2) we obtain d x F = 0, d p F = 2 and d q F = — 1, so that ao = bo = 0, 
Co = 1 and do = 2, therefore from Proposition 1 the local approximated solu- 



tion of first order is y(x) 



c , which is the exact local (global) solution of the 



original equation, satisfying the condition y(0) = 1. 

If F € C 1 ([a, 6] x R 2 ), F(xo,Po,qo) = and d p F(xo,Po,qo) 0, in a similar 
fashion the first order differential equation F(x, y(x), y'(x)) has the following 
local approximation of first order. 



y{x) = Po 



d x F 
d p F 
d x F 
d P F 



(x a ,p ,qo) 



(x ,po,q ) 



x 



d q F 

-g-j;{x ,Po,qo) 
d F 

-^(xo,po,go) 



qo - 
y\x) 



(25) 



If d q F(xo,Po, qo) = 0, the local approximated solution is simply 



3/0*0 = Po 



d x F 
dnF 



x 



d x F 
d p F" 



(26) 



while in the case d q F(xo,po, qo) 7^ one can verify that the local solution is the 
same as (|21f) . So the following statement holds: 

Proposition 2 In the hypothesis of Theorem 2, the ODE F(x,y(x),y'(x)) = 
has a local solution, approximated to first order, of the form 

y(x) = a Q + b x if d q F(x ,Po, qo) = (27) 

y{x) = a + box + coe doX if d q F(x ,Po, qo) £ (28) 

where ao, bo, cq and do are constants depending from (a;o,Po>9o)- 

Example 2 bis. For the ODE y'(x) 2 +y(x) 2 — 1 = 0, being F(x,p, q) = p 2 + q 2 — l, 
at point (5,1,0) we have d x F = 0, d p F = 2 and d q F = 0, therefore a = 1, 
bo = 0, Co = and do = 2. From previous proposition, a local approximated 
solution of first order is y(x) — 1. This result, compared with the exact solution 
y(x) — sin(x), is not satisfactory. A greater order of approximation is required. 

If F £ C n ([a, b] x K 2 ) for a suitable integer n > 0, F(xo,Po, qo) = and, 
e.g., d p F(xo,po,qo) ^ 0, it is possible consider higher order derivatives for the 
implicit function ip, as (II 81) or 



d qq ip(x ,qo) = - 



(dgip) d pp F + 2d q i> d 2 F + d 2 F 



d p F 



(29) 



G 



Then, from Taylor expansion of ip, we can write the following approximation of 
order (1,2) (1 for re-variable, 2 for g-variable) for an ODE: 



y(x) =p + 
1 



2d p F 



d p F 



2d x Fx + 2d q Fq 



d x Fx + d q Fq 



'[d q F] 2 d 2 F B q F 



[dp Ff 



d p F pq 



dl q F + d 2 qq F\ql 



+ 2^d pq F + d qq F \q \y (x) 



[dp FY 



d p F pq 



2d p F 



dp q F + d qq F I [y'(x)Y 



[dpF] »p* 
Expressions like this gives more accurate local solutions. 



Example 2 ter. For the differential equation of Example 2 we have d x F(^, 1, 0) = 
0,F(f ,1,0) = d 2 pq F(%,l,0) = 0, 5 p F(f ,1,0) =a p 2 p F(f ,1,0) = ^F(f,l,0) = 
2, so from previous formula the approximation is 

y{x) = l-\[y\x)f (30) 

The solution to this non linear equation, with condition y(-|) = 1, is (computing 
it using DSolve[{y[x]==l-l/2*y' [x] ~2,y[Pi/2]==l},y[x] ,x] of Mathemat- 
ica) 

, . 7T 2 7T 1 9 . . 

u(x) = 1 1 — x x (31) 

V ; 8 2 2 V ' 

This expression is the second order truncated Taylor expansion of sin(x) cen- 
tered at x = ? , so we have obtained a good approximation for the local solution 
found in Example 2. Note that, being in this case d x F = d^ x F — 0, from 
<(T7|) and l(T%|) follows that d% x ip{^,0) = 0, hence second order derivative respect 
a:-variable is not useful for a more accurate expression of the local solution. 



4 Relationship with series expansion of solution 

The last example seems establish a relation between symbolic expression of 
local solution found using implicit function technique and powers expression of 
solution found using series expansion of the unknown function in the original 
ODE. Let us consider the following symbolic expansion of second order, centered 
at xq = J, for the first member of the original equation y(x) 2 + y'{x) 2 — 1 = 
(using the Mathematica command ode = Series [y [x] ~2 + (D [y [x] ,x] ) ~2 - 
1, {x, Pi/2, 2}], see 0): 
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,7!\o //7T 



v3 2 + 2/(77) -y'U) vG) + v"0 (^-^) 



•2' a v 2' " v 2 



,7I\ „,7I\ 



'2' " v 2' 



»'(ff + /(f)(»(f) + /(f)) +9 '(fy»(f)](-| + 

Then we solve the algebraic system formed by the expanded differential equation 
and the two conditions = 1, ?/'(§) = (command Solve [{ode == 0, 

y[Pi/2] == 1, (D[y[x],x] /. x->Pi/2) == 0}]). The result consists of 
two possible choices: j/"(f ) = and £/"(§) = — 1. The expansion of the solution 
is y(x) = y(f ) + y'(Z)(x - f ) + |y"(f )(x - f ) 2 ; for y"(f ) = we obtain the 

trivial solution y = 1, for ^"(f) = —1 the solution is y(x) — 1 — ^ + t^x — ^x 2 . 
Therefore the series expansion method is computationally equivalent to the 
method of implicit function, but not logically: the first is based on a direct 
expansion of the unknown function in the original differential equation, while 
the second is based on an approximation of this ODE for finding a symbolic 
solution of the new equation. 

5 Finer examples 

Example 3. Let us consider the non linear ODE (see 

y'(x) = -3sin(x)[y{x)]* (32) 
This equation is integrable by elementary functions: 

27 

v{x) = ~rr^ — (33) 

[C + 6COS [X )\ 

The associated function F is F(x,p, q) — —3sin(x)p3 —q, and consider the point 
T = (f , 1, -3^) where F is null. The following identities hold: 

d x F = -3cos(x)p3, d p F = -4sin{x)pi , d q F(T ) = -1 (34) 
d x F(T ) = -^, d p F(T ) = -2y/3, d q F(T Q ) = -1 (35) 

We want a Taylor expansion for the implicit function q = <p(x,p) of order (3, 1), 
that is of third order respect to x-variable and first order respect to p-variable. 
We use expressions Ijl7|l and l|18|l for 0, and the following formula for the third 
derivative: 



3d x ^d 2 xx c^d 2 p F + [d x 4>? d 3 ppp F + 3d 2 xx cbd 2 xp F 



d q F 



(36) 



3 [d x <t>f dl pp F + 3d x ^ xp F + d 3 xxx F 



d q F 
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Computation gives the following values for the derivatives at T : 



d x <f> 



d 2 



3 

"2' ~ XXT ' 2 
so the new approximated ODE is 



,V3 



-2-\/3 (37) 



2 2( a; -|) + 4^(^-I) 



(38) 



Solving this equation with conditions &/(■§) = 1, J/'(§) = — 3-^2 (e.g with com- 
mand DSolve) one obtains a solution of the form y(x) = e ax (b + cx + dx 2 
which Taylor expansion of third order centered at -| is 



ex" 



(39) 



27 



y(x) = 1 + (tt - 3x) + ^ (tt - 3x) 2 + (tt - 3x) 3 

The exact solution of the original differential equation is y(x) = — — 
which Taylor expansion of third order centered at ^ is 

y(x) = 1 + ^ (tt - 3x) + ^ (tt - 3x) 2 + ^= (tt - 3x) 3 

Therefore the implicit function method gives a very good approximated sym- 
bolic solution. 



|+3cos(x)J 



(40) 




Figure 1: Comparison between exact and approximated solution in the interval 
[tj — 1, j + l] ; the exact solution is dashed. 
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Example 4- Consider a spherical bubble of a compressible gas immersed in an 
incompressible nonviscous fluid. If this fluid is flowing in a vane of a mechanical 
structure as rotating pump, presence of bubbles (cavitation) can give increment 
of noise and reduction of efficiency. Hence it is important to estimate the time t c 
necessary for the collapse of the bubble, assuming that gas pressure is neglegible 
compared to external fluid pressure pf. If y(t) is the radius of the bubble, it is 
possible to derive from Navier-Stokes equations the following ODE (see [H]) 

2 y{t)y"{t)+y'{tf = -\ P -j- (41) 

6 6 p 

It is convenient to express this ODE as a first order equation, which is equivalent 
if one doesn't consider constant solutions y(t) = k (see [§] for details): 

y(tfy'{t) 2 = \ P -l[Rl-y(tf] (42) 
o p 

where p is the density of the external fluid and Ro is the radius of the bubble 
at time t = 0. Without considering the trivial solution y(t) = Ro, which is 
not physically admissible, this equation is not symbolically integrable using 
elementary functions. 

For simplicity, but without lack of generality, suppose p = 1 and the pressures 
normalized with respect to the value of pf, so that pf = 1. Let then F(t,p, q) = 
|i?o — (| + q 2 ) P 3 the associated function F defining the implicit solution. From 
F(0, i?o,9o) = follows q = 0, as required from the physics of the problem, 
therefore the point (to,Po, <?o) is To = (0, Ro, 0). The following formulas hold for 
first and second derivatives of F: 



d t F = 0, d p F = -(2 + 3q 2 )p 2 , d q F = -2p 3 q, 
d 2 pp F = -(4 + 6q 2 )p, d 2 qq F = -2p 3 , d 2 q F = -6p 2 q 

from which we have d p F(T ) = -2i^, d q F(T ) = d% g F(T ) = 0, d 2 p F(T ) = 
-4i? and d 2 q F{T ) = -2i?g. Then, using (|16|) . i|17[l an d (|29|l . we can use the 
following approximation of order (1,2) (1 for y, 2 for y'), centered at (0,i?o,0), 
for g2J : 



y(t) = Ro 



(43) 



This equation is a non linear ODE, but it is integrable by elementary functions. 
Its particular solution with initial condition y(t) = Rq is (eg using command 
DSolve[{y[t]==R_0(l-l/2 y' [t] ~2) ,y [0] ==R_0} ,y [t] ,t]) 



Rn 



from which we have our estimate for the time t c of collapse: 

t c = V2R Q 



(44) 



(45) 
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From this formula one can observe the importance of the initial radius of the 
bubble and consequently develop suitable actions. It is interesting to compare 
this result with the estimate from a numerical resolution of the original ode l|41|l 
(NDSolve [{2/3y [t] y " [t] +y ' [t] ~2==-2/3 , y [0] ==R_0} , y [t] , {t , , . 12}] ). 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 



Figure 2: Comparison between symbolic approximated solution of the reduce ODE 
(continuous line) and numerical solution of the original ODE (dashed line); case pf 
. f> I Rq — 0.1; time is in horizontal axis, radius in vertical axis. 



The figure 2 shows that symbolic solution of the reduce ODE is a good approx- 
imation of the implicit solution of the original equation. 

6 Technical notes 

We have use internal commands DSolve and NDSolve of Mathematica ver. 5.2 
for symbolic and numerical resolution of ODEs cited in the examples. Intc- 
grability of these equations using elementary functions has been tested using 
DSolve (see [5]). 
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